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' Abstract 

in ! 

Q\ . A numerical search for bosonic superconducting static vortex rings in a 

U (I) a xU(l)w model is presented. The fate of these rings without current, 
is to shrink due to their tension until extinction. The superconductivity 

\Q . of the loop does not seem to prevent shrinking. Current quenching takes 

place before stabilization. 

a: 

Qh| 1 Introduction 

In a series of papers [U [2| classically stable, metastable quasi-topological do- 
main walls and strings in simple topologically trivial models, as well as in the 
two-higgs Standard Model (2HSM) were studied. They are local minima of the 
energy functional and can quantum mechanically tunnel to the vacuum, not be- 
ing protected by an absolutely conserved quantum number. In [3j a search for 
spherically symmetric particle like solitons in the 2HSM with a simplified higgs 
potential was performed without success. Although the existence of spherically 
symmetric particle-like solitons in the 2HSM has not been ruled out, we shall 
here look, instead, for axially symmetric solutions in a similar system. 

Consider a model with superconducting strings [2l Sj. Take a piece of such 
string, close it to form a donut-shaped loop and let current in it. A magnetic 
field due to the supercurrent will be passing through the hole of the donut (figED). 
The energy of the loop has, a term proportional to the length of the string and 
will tend to shrink the radius of the donut to zero and the ring to extinction. 
However, another force opposes this tendency. Namely, as the loop shrinks, the 
magnetic field lines are squeezed in the hole, since, due to the Meissner effect, 
they cannot leave the loop. They are trapped inside the hole of the donut, oppose 
further shrinking and might even stabilize the string. 
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Figure 1: A x — y profile of the superconducting ring (left) as well as a y — z profile 
(right) where one can view how the mechanism we propose, against shrinking, works. 

This, as well as other arguments [7j are inspiring but not conclusive. The 
magnetic field will not be trapped inside the loop if the penetration depth of the 
superconductor is larger than the thickness of the ring. Also, once the magnetic 
field gets strong it can destroy superconductivity and penetrate [TTJ. Further, 
there is a maximum current a superconductor can support (current quenching). 
This sets a limit on the magnetic field one can have through the loop, and this 
may not be enough to stabilize it. Thus, the above approach may work at best in a 
certain region of the parameter space, depending also on the defect characteristics 
[2]. The purpose of this work is to apply the above straightforward idea to search 
for string loops [l3]-[22] in a U(l) x U(l) gauge model and to determine the 
parameter space, if any, for their existence and stability. Another interesting 
subject is to have a rotating ring. The rotation is another extra factor which 
could help the ring to stabilize. This work was done with success both analytically 
and numerically on [6] where vortons are exhibited. Another recent example of 
rotating superconducting electroweak strings can be found on [8], while for a 
review on electroweak strings, the reader should also check |9|. Finally, a work 
on static classical vortex rings in SU (2) non-abelian Yang-Mills-Higgs model can 
be found on [TO] . 

2 The model 

The Lagrangian describing our system is: 

£ = -\*%> - \ w l + lA^I 2 + \d^\ 2 - u{\4>\, H) (l) 
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where the covariant derivatives are D^ip = d^ip + ieA^ip, D^ifi = <9 M + iqW^, 
the strength of the fields are F^ v = d^A v — d u A^, = dpW v — d-JW^ while e 
and q stand as the relevant U(l) charges. We choose the potential U 

U(\4>\, M) = |(|0| 2 - vlf + f (M 2 - vl) 2 + f |0| V| 2 - f, 2 4 (2) 

The vacuum of this theory is |0| —Vi^O, \if)\ = and breaks U(l)w x U(1)a — ► 
U(1)a, giving non-zero mass to W. The photon field stays massless. There, 
U(vi,0) = 0. The vacuum manifold M. in this theory is a circle S 1 and the first 
homotopy group of M. is tti(A4) = vri(S' 1 ) = Z which signals the existence of 
strings. In regions where \<p\ = 0, the field |^| is arranged to be non-vanishing 
and U(l)w x U(1)a — > U(l)w- Thus, U(1)a — > 1 and one may generate an 
electric current flowing along regions with vanishing \<f>\. Hence, this theory has 
superconducting strings [4j. The vacuum of the theory leaves unbroken the elec- 
tromagnetic U(1)a- For g 3 vf > g 2 v 2 2 this vacuum is stable, while g\v\ > g 2 v\ 
ensures that it is the global minimum of the potential. The mass spectrum is 

m A = 0, m w = qv 1 , m 2 4> = g 1 v\, m\ = - (g 3 v 2 - g 2 vf) (3) 

3 The U{1) A x U(l)w model: Search of superconducting vortex rings 

Configurations with torus-like shape, representing a piece of a U(l)w —> 1 
Nielsen-Olesen string [TJ], closed to form a loop, are of interest in this search. 
Thus, we will require <p to vanish on a circle of radius a (the torus radius) 
<f)(p = a,z = 0) = 0. At infinity (p — ► oo, z — > oo), we have the vacuum of 
the theory. This translates to \(j)\ — ► Vi, — > 0. The ansatz for the fields is: 

<j>(p,<p,z) 
ip(p,<f,z) 

A{p,(p,z) 
W(p,<f,z) 

where M, N are the winding numbers of the relevant fields and p, </?, z are the 
cylindrical unit vectors. We use cylindrical coordinates (t,p,<p,z), with space- 
time metric g^ v = diag(l, — 1, — p 2 , — 1). We define Q(p, z) = arctan(j;/(p — a)). 
We work in the A = = W° gauge. Especially for the gauge fields, we suppose 
the above form based on the following reasonable thoughts: The W-field is the 
one related to the formation of the string thus, it exists in the constant-y? plane. 
This means that in general its non- vanishing components will be W p and W z . The 
A-field is the one produced by the supercurrent flowing inside the toroidal object. 
The current is in the (p direction thus, we in general expect the non-vanishing 



= F(p,z)e iM@ip ' z) 

= <p 

p 

= WJp,z) p + W z (p,z) z 



3 



1-1 



Figure 2: A wider view of the "Mexican hat" potential (eq|5]) for (gi,g2,u) = (10,8,1). 

component to be A v . Finally, as it concerns the amplitude of all the fields, we 
expect that it is independent of ip due to the axial symmetry of torus. 

With the above ansatz, the energy functional to be minimized takes the 
form: 



E = 2rt Vl [ pdp [ dz ((cyy 2 + (d z A^) 2 ) + 

JO J-oo V J 



+ (9 p P) 2 + (c) z P) 2 + ^(e^ + iV) 2 H 



+ (d p Ff + (d z F) 2 + -(d p W z - d z W p f + 

+ {[qW p + Md p Q) 2 + {qW z + M^9) 2 ) F 2 + U{F, P) 



(4) 



and the potential U (figE|) can be written as follows: 



U(F,P) = y i{F 2 -l) 2 + 9 -l{P 2 -u 2 ) 2 + 9 -^F 2 P 2 - |« 4 



(5) 



where u = Vijv\. This is the energy functional we use for our computations. The 
conditions to be satisfied by the parameters become: 



9i > 92U 4 , & > £ 2 w 2 



(6) 



The gauge fields A, W have magnetic fields of the following form: 



V x A = B , 



1 / dAu, . dA u 



p\ dp 



z — 



dz 



V x W = B 



w 



dW z _ dW p 
dp dz 
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The field equations follow: 



d 2 p F + d 2 z F+^-^(F 2 -l)F 



93p2 F 



, zcos 2 Q \ ( „r „, cos 2 © X 
qW p - M— — I + I qW z + M 



(p - af 



(p-a) 



F = 



+ 9 2p + ^_ _ ( ^ + N } 2 P _ 92 ( p2 _ ^ p _ 9_l F 2p fJ 



#X - d z 8 p W z - 2qF 2 \qW p -M 



zcos 2 Q 

Ip^ 2 

d 2 p W z + - [d p W z - d z W p ) - d p d z W p - 2qF 2 I qW z + Mj^- 
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We can also write down the currents associated with field, namely j% and 
jf and the total current X^ out of these as well as the supercurrent X^ associated 
with the if; field. These are 



\l{jf) 2 + u 



<t>\2 



T 



2eP 2 



(eA v + N) 



(7) 



where 
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-2qF 2 (qW p + Md p B), jf = -2qF 2 (qW z + Md z Q) 



(8) 



As explained in the introduction, Meissner effect is of crucial importance 
for the stability of the torus-like object. The magnetic fields produced by the 
supercurrent X^ penetrate into the toroidal defect in a distance dictated by the 
penetration depth. In this theory, the U{l)x symmetry breaks inside the string 
and the photon acquires mass m 2 A = e 2 < P > 2 where < P > is the expectation 
value of the charge condensate in the vicinity of the string core. There is the 
superconducting sector of the defect. The penetration depth is A = ^- = e< p > . 
But < P >= P m ax < u thus, one can have an estimate for A, with a lower bound 
for its value. This is A > — where the equality holds when P ma x = u. On the 
other hand, one can also compute an upper bound for the thickness of the defect. 
We know that ta, = — = -t=. If our concern is to search for stable rings, a 
reasonable step is to demand the penetration depth to be smaller than the string 
thickness which means 

1 1 



A < ?> 



eP 



< 



e P max > 91 



(9) 



This is the condition needed in this case. From the above condition, we get the 
diagram shown in figj3] where one can see the area where it's more possible to find 
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Figure 3: For different values of ^/g%, we plot the lower bound over which the condition 
(JH) is satisfied. The plot is P m ax vs. e and we plot the function P max = -^p-- 



stable solutions if they exist. The numerical results we present later, are what 
we found while searching inside this region. 

A way to derive virial relations is through Derrick's theorem. The virial 
relation for the field of our model, must have the constraint <fi K (p — a,z) = 
0(p = a, kz) = 0. Consider the rescalings p — > p, z — > kz, F k — > F, P K — > P, 
Atp K — > A^,, Wp )2;K — > kW P)Z . Then, we find the minimum through the relation 
W- = when k — 1. The virial relations we use in our model, in order to check 

OK ' 

our numerical results follow. If we define 



/oo poo 
pdp I dz 



~{d p W z - d z W p f + ^(d z A v ) 2 + {d z P) 2 + {d z Ff + 



+ 2F 2 qW p (qW p + Md p Q) + (qW z + Md z Bf 



I2 = —2nvi j pdp I dz 



^- 2 {d P A v ? + {d p Ff + {d p Pf + ^{eA v + Nf + 
zp p 



+ {d z W p ){d p W z - d z W p ) + F 2 {qW p + M<9 P 6) 2 + (qW z + Md z Q) 2 + 



+ ^( J p2_ 1 )2 + ^ ( p2_ M 2 ) 2 + ^3 F 2p2_^4 



we must have h + h = 0. We define the index V = ijnj^r and we want its value 
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to be as small as possible. We can derive many other virial relations by assuming 
generally for a field 0, the "double" rescaling <f>(x) — > K(f>(p,x) and then demanding 
|^| re=1=jU = = |p| re=1=M . For example, consider the following rescalings p — > p, 
z -> fiz, F K -»• P, P K -> f«P, A Vk -> VT P)2k -> W p ,«. We have 



J 3 = 27rt>i y pdp / dz 



2p 



2 (d z A^) 2 + (d z P) 2 + {d z Ff + 



+ 2P 2 ( M^6( 9 W a + M<9 2 6) ] + -(d p W z - <9 2 W P ) 



I4 = —2-kvx I pdp I dz 



^(<Vy 2 + («9 p P) 2 + (d p F) 2 + ^(e^ + N) 2 + 
zp p 



+ d p W z {d p W z - d z W p ) + P 2 {qW p + M<9 P 9) 2 + (gW* + Md z Qf + 



4 v ; 4 v ; 2 4 



together with 



POO 

/ pdp 

/O j 


POO 

/ 

' — 00 


POO 

/ pdp 

/0 j 


POO 

/ dz 

' —00 



2P 2 

2(«9 p P) 2 + 2(d z P) 2 + —(eA v + N) 2 

P 



^2 D 2 



where, as above, we must have ^3 + I4 = = 1$ + I§. 



3.1 Numerical results 



Energy minimization algorithm is used to minimize the energy functional of (j4j). 
The algorithm is written in C. One can find details about the algorithm on page 
425 of [5] but, briefly, the basic idea is this: Given an appropriate initial guess, 
there are several corrections to it, having as a criterion the minimization of the 
energy in every step. When the corrections on the value of the energy are smaller 
than ~ 10~ 8 the program stops and we get the final results. A 90x20 grid for 
every of the five functions is used, that is, 90 points on p-axis and 20 on z. Here, 
we begin with fixed torus radius a. Then, the configuration with minimum energy 
for this a is found. Other values of a are chosen as well and the same process 
goes on until we plot the energy vs. the torus radius E(a). It would be very 
interesting to find a non-trivial minimum of the energy (on a m i n 7^ 0), which 
would correspond to stable toroidal defects with radius a m i n . 
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The initial guess (figSl) we use for our computation is: 



F(p,z) = tanh((p - a) 2 + z 2 ) M l 2 

P{p,z) = tanh(p 7V )(l-tanh((p-a) 2 + z 2 ] 



Afpip^z) = tanh 

W p {p,z 



2 



' ' ' ' ' e \((p - a) 2 + z 2 ) 2 

Mz cos 2 q( {p-a) 2 + z 2 



W z (p,z) 



q(p-a) 2 y(p-a) 2 + z 2 + (a 2 /A) 
M cos 2 ef (p-a) 2 + z 2 



q(p-a) \(p-a) 2 + z 2 + (a 2 /4) 

where £ a constant. This initial guess also satisfies the appropriate asymptotics 

• near p = 0: 

F^O, P~p N , A^p 2 j{z) (10) 

• near (p = a, z = 0): 

F~p M/2 , W p = = W z (11) 

• at infinity: 

F ~ l-0(e - ^), P~C(e-v / ^+ 15 ) 

W P — ^eU + ^e-^), W 2 ~-— 9 z eU + G(e-^) (12) 

q q 

where p = (p — a) 2 + z 2 . 

For fixed torus radius (i.e. here a ~ 2.2) we present a typical graph of the 
final configuration of the lowest energy (see figEl). We also present the plot of 
the energy of the system vs. the radius of the toroidal object which reveals the 
instability of the system as there is no non-trivial minimum (see figE]). 

From the results, we can point out a few things. Firstly, the greater the 
value of e we use, the stronger the supercurrent becomes. Secondly, the greater 
the value of e we use, the lower the radius a where the supercurrent quenches 
(figE]). These are expected as the increase of e makes the condition of eqj9] 
stronger, something which means that the mass of the photon increases and 
the penetration depth decreases at the same time. It is also reasonable that 
a stronger current can "defend" the defect, against the magnetic field, a little 
longer. Another observation is that as the difference g± — g>2 decreases, that is to 
say, as g2 increases and becomes close to g±, we see that the maximum current 
increases and the quenching takes place again at lower a (figE]) . This is expected 
as one can see from the potential (eq]5]) of the energy functional, because the 
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/\(p'0) 

1 2 3 4 5 6 

P 

Figure 4: A typical plot of the initial guess we use for the five fields for the lowest 
winding state M = 1, N = 1 on the z = plane. 

stronger the coupling g 2 is, the more important the relevant term g 2 (P 2 — u 2 ) 2 /A 
becomes. The latter has as a consequence, the increase of P max which counteracts 
the effects from the increasing g 2 coupling. 

The parameter space where we searched, starts from gi — 4 (m^ = 2). In 
order to search the model, we reached values around g\ = 30 (m^ « 5.48) over 
which, e has to be very large in order to respect the condition ([9]). There is also 
the fact that great values of e in general is an unwanted feature since we use 
semiclassical approach. As it concerns the other couplings we have g% = g\ and 
g 2 = gi - k with 0.5 < k < 8, (u = 1). 

4 Explanation concerning the instability of the vortex ring 

From the condition in eqJH], it is understood that we are enforced to lower gi 
as much as possible and/or increase e. But these steps are not as easy as they 
might look. There are some limits. The lower bound on the value of g\ has a 
reasonable explanation. For low values of g\, the coupling g 2 is also low (because 
9i > 92U 4 ). Now, when g 2 is small enough, the changes on the term g 2 (P 2 —u 2 ) 2 /4: 
are unimportant for the energy, in comparison to the term (d p P) 2 . In that case, 
the lowering of the last term minimizes the energy, something which means that 
P — » 0. Indeed, this is computationally observed. There is also a lower limit on e 
which is reasonable because the lowering of e results to an increasing penetration 
depth. 

We searched for stable objects for values over these limits described above. 
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P 

Figure 5: Typical plot of the final configuration of fields on the z = plane. Parameters 
are M = 1, N = 1, e = 5, q = m w = 2, gi = 14, = Jgi = 3.74, g 2 = 12, g 3 = 14, 
u = 1, rruj, = a/ (53 — g2U 2 )/2 = 1, vi = 7.5 • 10~ 3 . Energy E = 0.72 and virial is 
1.5 • 1(T 3 . 

The computational results are exhibited in figsJMHland as we see, these objects are 
unstable. We argue that the explanation for the instability is current quenching 
and that, only for high values of e. For values of e, of the order of 1, we have, 
according to eqEJ that the penetration depth is much bigger than the string 
thickness thus, stability is out of the question. The latter is also computationally 
observed. 

Now, we base our aspect about quenching on qualitative as well as quantitative 
arguments. We observe that as the torus shrinks, the supercurrent suddenly drops 
to zero which signals the destruction of the defect. Just before the sharp drop, we 
notice that the supercurrent rises with increasing rate. This must be due to the 
resistance the torus meets from the magnetic lines as it shrinks. One can observe 
that as the supercurrent increases and the condition of eqj9] is satisfied at the 
same time (i.e. see dashed and dotted line of figEl), suddendly the current is lost. 
This can be explained only through current quenching. The above phenomenon is 
not observed when eqj9]is not satisfied (i.e. see dotted line of figlZl). There, as the 
magnetic lines penetrate the ring, they meet almost no resistance since A is much 
greater than r^. Another observation which supports our quenching argument is 
that, rough estimations on the maximum current a string can sustain, lead us to 
the following formula (see p. 129 of [11]) which makes a small overestimation in 
order to have an upper limit: 

TL. < vW (13) 
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Figure 6: The top graph is the energy vs. the torus radius E(a). The middle graph is 
the supercurrent vs. a. On that graph one can clearly see current quenching. The 
bottom graph is the quantity e 2 P% iax vs - a ( or m A vs. a) where one can see the area 
in which the condition of eq|9] holds (lines over the gi-limit line). The resistance from 
the magnetic field can be seen as a sharp increase on the supercurrent. Dotted lines 
are for e = 6, dashed for e = 8 while dashed and dotted for e = 10. All plots are for 
(M, N, u, gi, 52, $3)=(1, 1, 1, 14, 12.5, 14). 



where a = J J dpdz P 2 . In figsEHZl the maximum value of the supercurrent 
is close to the limit of the estimation of eqfT3l The table below gathers the esti- 
mated Tf st (according to eq JT3l) and the computed maximum supercurrent (Zf om ) 
for the parameters of these figures. 



e 


9\ 


92 


#3 


7-0 

M est. 


com. 


6 


14 


12.5 


14 


4.6 


4.3 


6 


14 


13 


14 


5.2 


5.0 


8 


14 


12.5 


14 


6.7 


6.4 


10 


14 


12.5 


14 


8.0 


7.2 



Finally, one can make an estimation of the value of the supercurrent which 
would stabilize the ring, namely %f tab . This can be done as follows. As ex- 
plained in the introduction, there is the tension of the string which shrinks 
the loop and the magnetic field which opposes this tendency. When the ring 
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Figure 7: The top graph is the energy vs. the torus radius E(a). The middle graph is 
the supercurrent vs. a. On that graph one can clearly see current quenching. The 
bottom graph is the quantity e 2 P% iax vs - a ( or m A vs. a) where one can see the area 
in which the condition of eq|9] holds (lines over the gi-limit line). The resistance from 
the magnetic field can be seen as a sharp increase on the supercurrent. Dotted lines are 
for gi = 12, dashed for <?2 = 12.5 while dashed and dotted for 52 = 13. All plots are for 
(M,N,u,e, 9l ,g 3 )=(l,lA, 6,14,14). 



is stabilized we have E tension = E magnetic . Here, E tension ~ 1 and E magnetic = 
2irvi J pdp J dz(B\/2). Without any calculation, one can point out that, since 
the total energy in the quenching radius is below E tension = 1, then the E magnetic 
which is a fraction of it, would be even smaller. Recall that Ba oc X^, which means 
that we need a lf tab which is well above the maximum current we can have inside 
the defect. Calculations of the magnetic energy are in agreement with the above 
observation and place its value around E magnetic ~ 5 • 1CT 3 << E tension ~ 1. This 
translates to the following conclusion: 2^ tab > 10 ■ X^ st . 

Thus, we find out that the current needed for stabilization, is at least ten 
times bigger than the value of the maximum current our string can sustain. We 
also observe that our computational maximum current values are close to the 
theoretical estimations about quenching. This means, that we will have current 
quenching as an "obstacle" towards stabilization, since the supercurrent will not 
be enough in order to create the magnetic field needed. 
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Figure 8: The plot of the energy of the system vs. the radius of the torus for four differ- 
ent sets of parameters. Dotted is for (gi, g2, <?3)=(14, 12.5, 14), dashed is for (18, 15, 18), 
dashed and dotted is for (25,20,25) while dashed with two dots is for (30,24,30). For 
all data sets we have (e, M, N, u)=(10, 1, 1, 1). As one can observe, there exists no 
minimum. 

5 Conclusions 

Future experiments on LHC could answer whether metastable particle-like soli- 
tons exist in MSSM or 2HSM or not. On [3] there is a search for spherically 
symmetric solitons in the frame of the 2HSM but with a simplified potential. 
Here we search for axially symmetric solitons which, if stable, will have a mass of 
few TeV [2|. We considered the U(1)a x U(1)w model, where the existence of the 
vortex is ensured, for topological reasons. There, we searched for stable toroidal 
strings. We present and analyse our observations. This paper tries to answer to 
expectations having to do with observations of stable axially symmetric solitons 
which would be possible to detect in later experiments of LHC. For relatively 
small values of e (~ 1), the system seems to have no stable vortex rings. In fact, 
this instability is present even in other parameter areas where we searched (i.e. 
e > 6 see figs J6|7l) . We explain why we believe that the main reason of instability 
is current quenching. 
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